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Abstract 



We analyze the effects of charged lepton corrections and renormalization group (RG) running 
on the low energy predictions of theories which accurately predict tri-bimaximal neutrino 
mixing at the high energy scale. In particular we focus on GUT inspired see-saw models 
with accurate tri-bimaximal neutrino mixing at the GUT scale, in which the charged lepton 
corrections are Cabibbo-like and give rise to sum rules valid at the GUT scale. We study 
numerically the RG corrections to a variety of such neutrino mixing sum rules in order to 
assess their accuracy and reliability when comparing them to future low energy neutrino 
oscillation experiments. Our results indicate that the RG corrections to neutrino mixing 
sum rules are typically small (less than one degree), at least in the examples studied with 
hierarchical neutrinos. 

1 Introduction 

Perhaps the greatest advance in particle physics over the past decade has been the discovery of 
neutrino mass and mixing involving two large mixing angles commonly known as the atmospheric 
angle 023 and the solar angle 6i2- The latest data from neutrino oscillation experiments is 
consistent with the so called tri-bimaximal (TB) mixing pattern [Tj , 



where PMaj is the diagonal phase matrix involving the two observable Majorana phases, and 
there were many attempts to reproduce this as a theoretical prediction [2l[3llll[5l[6l[71[8l|9l [10] . 
Since the forthcoming neutrino experiments will be sensitive to small deviations from TB mixing, 
it is important to study the theoretical uncertainty in such TB mixing predictions. 
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The question of how to achieve TM mixing has been the subject of intense theoretical 
speculation. In theoretical models one attempts to construct the neutrino and charged lepton 
mass matrices in some particular basis. There are two particular bases that have been used 
widely in the literature for this purpose, as follows. The first basis is the flavour basis in which 
the charged lepton mass matrix is diagonal, while the neutrino mass matrix takes a particular 
form such that is results in TB mixing. The second basis is a particular basis first introduced 
by Cabibbo and Wolfenstein in which both the neutrino and charged lepton mass matrices are 
non-diagonal, but in which the charged lepton mass matrix is diagonalised by a "democratic 
unitary matrix" involving elements of equal magnitude but differing by a phase u) = 27r/3. Such 
a Cabibbo- Wolfenstein basis is particularly well suited to models of TB mixing based on the 
discrete group However in other classes of models, one attempts to work in the flavour 

basis and to derive TB mixing purely from the neutrino sector with the charged lepton matrix 
being diagonal, for example using constrained sequential dominance (CSD) [2]. 

However, when attempting to derive TB mixing in the flavour basis in realistic models arising 
from Grand Unified Theories (GUTs), it is observed that, although TB mixing may be accurately 
achieved from neutrino mixing, in practice the flavour basis is never accurately achieved, i.e. 
the charged lepton mass matrix is never accurately diagonal. Instead, in such GUT models, 
the charged lepton mass matrix often resembles the down quark mass matrix, and involves an 
additional Cabibbo-like rotation in order to diagonalize it. In such models, then, TB mixing 
arises in the neutrino sector, but with charged lepton correction giving deviations [12 . It turns 
out that such Cabibbo-like charged lepton corrections lead to well defined corrections to TB 
mixing which can be cast in the form of sum rules expressed in terms of the measurable PMNS 
parameters. One example is the neutrino mixing sum rule 612 — ^13 cos((5) « ^12 111 ESI E] j where 
5 is the observable Dirac CP phase in the standard parameterisation. Since such sum rules may 
be tested in future high precision long baseline neutrino experiments, it is of interest to know 
with what precision they are expected to hold theoretically. 

Such sum rules as discussed above are strictly only expected to apply at some high energy 
scale, whereas the neutrino experiments are performed at low energy scales. In order to compare 
the predictions of such sum rules to experiment one must therefore perform a renormalisation 
group (RG) running from the high energy (e.g. the Grand Unified Theory (GUT)) scale where 
the theory is defined to the electroweak scale Mz- RG corrections arise mainly from the large 
tau lepton and third family neutrino Yukawa couplings, and this leads to large wavefunction 
corrections in the framework of supersymmetric models. The running of neutrino masses and 
lepton mixing angles is very important and has been studied extensively in the literature |15) . 
In |16| fTT! a mathematica package REAP (http://www.ph.tum.de/ rge/) which solves RGEs 
and provides numerical values for the neutrino mass matrix and mixing angles was developed. 

In this paper we provide a first numerical study of the deviations from TB mixing due to 
the effects of both charged lepton corrections and RG running. We focus on GUT inspired 
models in which the charged lepton corrections are Cabibbo-like and in this case they may be 
cast in terms of sum rules valid at the GUT scale. In practice, then, we are interested in the RG 
corrections to these sum rules which may subsequently be reliably compared to experiment at 
low energies. We shall study a variety of neutrino mixing sum rules (arising from the deviations 
from exact tri-bimaximal neutrino mixing due to Cabibbo-like charged lepton corrections) and 
comment on their accuracy and reliability when comparing them to future low energy neutrino 
oscillation experiments. Most of the specific numerical results are inspired by a particular class 
of GUT-flavour models, namely the models in [H [B], which are precisely the type of models 
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in which the sum rules emerge in the first place. In the cases studied we find rather small 
corrections. For example the sum rule 9i2 — 0i3cos((5) ~ 35.3° becomes renormalized by about 
0.4° for large tan f3 = 50. Although most of the numerical results are based on a particular GUT 
motivated model, we also analyze a completely different type of model and find qualitatively 
similar results. This suggests that our results will apply to more general models based on the 
Minimal Supersymmetric Standard Model, extended to include the see-saw mechanism, with 
hierarchical neutrino masses. 

We emphasize that in this paper we only consider deviations from TB lepton mixing due to 
the combination of charged lepton corrections and RG running, and that in general there will be 
other sources of deviations that we do not consider. For example, exact CSD [2| itself will lead 
to some deviations since it does not predict precisely TB neutrino mixing due to corrections of 
order 777-2 / n^-s [H] • Another example of deviations to TB mixing are the canonical normalization 
effects discussed in [TH]. Clearly such additional corrections are quite model dependent and a 
phenomenological study of corrections to the TB mixing in the neutrino sector, in the flavour 
basis, has been made in [20]. In this paper we shall simply assume precise TB neutrino mixing 
at the GUT scale, and investigate the deviations due to charged lepton corrections and RG 
running only, ignoring other possible model dependent corrections. The RG corrections to TB 
mixing (but not charged lepton corrections or the resulting sum rules) were also considered in 

m- 

The layout of the remainder of the paper is as follows. In section 2 we give our conventions, 
including the TB deviation parameters that we will use. In section 3 we discuss the see-saw 
mechanism and show how TB neutrino mixing can be achieved. In section 4 we show how 
Cabibbo-like charged lepton corrections leads to neutrino mixing sum rules. We also present 
a numerical model of TB neutrino mixing adapted from a well motivated example of light 
sequential dominance used in the GUT flavour models of [21 [6] that we shall use in most of the 
remainder of the paper. We also show how sensitive these results are to non-Cabibbo like charged 
lepton corrections. In section 5 we numerically study the RG corrections to the various neutrino 
mixing sum rules which embody the charged lepton corrections to TB neutrino mixing. In section 
6 we explore a second type of numerical model adapted from heavy sequential dominance |18j 
and show that the results are qualitatively similar to the case of the first numerical model. 
Section 7 concludes the paper. 

2 Conventions 

2.1 The PMNS matrix in the standard parametrization 

The mixing matrix in the lepton sector, the PMNS matrix Upmns^ is defined as the matrix 
appearing in the electroweak coupling to the W bosons expressed in terms of lepton mass 
eigenstates. The Lagrangian is given in terms of mass matrices of charged leptons Me and 
neutrinos mjy as, 

C = -eiMeCR - -VLruLLvi + H.c, (2.1) 
The change in basis from flavour to eigenbasis is performed by, 

Ve^MeV^ = diag{me,mfj,,mr), Vy^niLLVj = diag{mi,m2,m3), (2.2) 
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The PMNS matrix is then given by, 

UpMNS = 

In the standard PDG parametrization, the PMNS matrix can be written as, 



(2.3) 



(C12C13 S12C13 siae"'" \ 

-C23S12 - Sl3S23Cl2e"' C23C12 - Sl3'S23Sl2e'^ S23C13 Phlaj, (2.4) 
S23S12 - Sl3C23Cl2e*'^ -S23C12 - Sl3C23Sl2e*'^ C23C13 / 

where S is the Dirac CP violating phase, and PMaj = diag{e^^ ,e^^ ,0) contains the Majorana 
phases ai, 02- The latest experimental values and errors for the three neutrino oscillation 
parameters are summarised inTablejl] \22\ l23l [M] . 



Parameter 


Best fit (°) 


2a n 


3a n 


O12 


34.44 


31.94- 37.46 


30.65- 39.23 


G23 


45 


38.05 - 52.53 


35.66 - 54.93 


ei3 


4.79 


< 10.46 


< 12.92 



Table 1: Best fit values, 2 a and 3 a intervals for the three- flavour neutrino oscillation param- 
eters from global data including accelerator (K2K and MINOS) and solar, atmospheric, reactor 
(Kam LAND and CHOOZ) experiments [22]. 



2.2 A parametrization in terms of deviations 

Another parametrisation of the lepton mixing matrix can be achieved by taking an expansion 
about the tri-bimaximal matrix. Three small parameters r, s and a are introduced to describe 
the deviations of the reactor, solar and atmospheric angles from their tri-bimaximal values |25], 



r 1 1 

■si3 = ^, si2 = ;^(l + s), S23=;^(l+a). (2.5) 

Global fits of the conventional mixing angles [24J can be translated into the ranges, 

< r < 0.22, -0.11 < s < 0.04, -0.12 < a < 0.13. (2.6) 

Considering an expansion of the lepton mixing matrix in powers of r, s, a about the tri- 
bimaximal form. One gets the following form for the mixing matrix to first order in r, s, a 



Umns 



1 



V6 



(1 - 



(1 -|- s — a -|- re* 



1 

V3 
1 , 



^ {l + s + a- re'^) 



-^(1-is + a + ire*^) 



Maj- 



(2.7) 
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2.3 Third row deviation parameters 



For later convenience, we also define the following parameters which express the deviation of 
the magnitude of the third row mixing matrix elements from their tri-bimaximal values: 



\{UpMNs)3l\ 
|(f^PAfAf5)32| 
\(UpMNs)33\ 



1 

1 
1 

71 



1 + ^3) 



(2i 



Hence from Eq 2.4 



6 
6 
6 



\/6 |S23S12 - Si3C23Cl2e* 
S23C12 - Sl3C23Sl2e 
\/2 |C23C13| - 1- 



1, 
- 1, 



I - S23C12 - Sl3C23Sl2e*'^| - 1, (2.9) 

We can also express the parameters in terms of the deviation parameters r, s, a from 

\l + s + a- re'^l - 1, 

\1- y + a+^re'^\-l, (2.10) 



Eqj2.7^s follows, 



6 
6 



1 



1. 



3 The See- Saw Mechanism and TB Neutrino Mixing 

3.1 The see-saw mechanism 

The see-saw mechanism provides an excellent explanation for the smallness of neutrino masses. 
Before discussing its simplest form, we first start by summarising the possible types of neutrino 
mass. 

One type is Majorana masses of the form rriLLVi^v'f^ where vl is a left-handed neutrino field 
and is the CP conjugate of a left- handed neutrino field, in other words a right handed 
antineutrino field [26j. Introducing right-handed neutrino fields, other neutrino mass terms are 
possible. There are additional Majorana masses of the form M^pii/^v'f^ and Dirac masses of the 
form ruLRUi^vji where i^r is a right-handed neutrino field and is its CP conjugate. 

The Majorana masses of the form m^^are strictly forbidden in the standard model, assuming 
only the higgs doublets are present. The reason for this is that heavy left-handed neutrinos would 
disturb the theory of weak interactions with W, Z bosons. For the simplest version of the see- 
saw mechanism, one can assume that the left-handed Majorana masses are zero at first, but are 
effectively generated after introducing the right handed neutrinos [27'. 

The right-handed neutrino does not take part in weak interactions with the W, Z bosons, 
and so its mass Mrr can be arbitrarily large. With these types of neutrino mass, the see-saw 
mass matrix is given as. 

In the approximation that Mrr ttilr [ Mrr may be orders of magnitude larger than 
the electroweak scale) , the matrix in Eqj3.1| can be diagonalised to give the effective left-handed 
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Majorana masses, 

friLL = mLRM~l^m[ji (3.2) 

These masses are naturally suppressed by the heavy scale Mrr. Taking niLji = My/ = 
SOGeF and Mrb. = Mgut = W^^GeV, we find tull ^ ld~^eV which is good for solar neutrinos. 
A right-handed neutrino with a mass below the GUT scale would be required for atmospheric 
neutrinos. 

The fundamental parameters which must be inputted into the see-saw mechanism are the 
Dirac mass matrix m/^/j and the heavy right-handed neutrino Majorana mass matrix Mjj/j. The 
output is the effective left-handed Majorana mass matrix m^^ as given by the see-saw formula 



in Eqj3.2|[28]. 

The see-saw mechanism discussed so far is the simplest version and it is sometimes called 
type I see-saw mechanism. In Pati-Salam models or grand unified theories based on 5*0(10), 
type I is generalised to type II see-saw, where an additional term for the light neutrinos is present 

EHl. 



3.2 Approximate TB neutrino mixing from CSD 

Sequential dominance (SD) is an elegant way of accounting for a neutrino mass hierarchy with 
large atmospheric and solar mixing angles |30| [3T| . Here we review how tri-bimaximal neutrino 
mixing can result from constrained sequential dominance (CSD) [2]. In SD, the atmospheric 
and solar neutrino mixing angles are obtained in terms of ratios of Yukawa couplings involv- 
ing the dominant and subdominant right-handed neutrinos, respectively. To understand how 
tri- bimaximal neutrino mixing could emerge from SD, we begin by writing the right- handed 
neutrino Majorana mass matrix Mrh in a diagonal basis as. 



Mrr « X . (3.3) 




Without loss of generality write the neutrino (Dirac) Yukawa matrix l^i? terms of the complex 
Yukawa couplings a,b,c,d,e,f,a',b',c' as 

(d a a' \ 
e h h' \. (3.4) 
fed) 

For simplicity we assume that d = 0. SD then corresponds to the right-handed neutrino of 
mass Y being the dominant term while the right- handed neutrino of mass X giving the leading 
sub-dominant contribution to the see-saw mechanism. 

where x,y £ a,b,c and x',y' £ a',b'c', and all Yukawa couplings are assumed to be complex. 
Light sequential dominance corresponds to 

Y <X <X'. (3.6) 
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Tri-bimaximal neutrino mixing, in which tan = 1; tan — ^ /v^ and ^J'g = corresponds 
to the choice, 

\a\ = \h\ = |c|, (3.7) 

\d\ = 0, (3.8) 

\e\ = I/I (3.9) 

e*b + f*c = 0. (3.10) 

This corresponds to constrained sequential dominance (CSD)[2]. Note that the analytic results 
for SD and CSD are accurate to leading order in 1712 /m^ [18:, so these conditions will not give 
rise precisely to TB neutrino mixing, and so in the numerical studies we shall need to perturb the 
CSD relations in order to achieve accurate TB neutrino mixing at the GUT scale, as discussed 
later. 



4 Charged lepton corrections 

4.1 Cabibbo-like corrections 
4.1.1 Sum rules 

In this paper we shall consider the case that TB mixing (Eq ]l.l[ ) applies quite accurately only to 
the neutrino mixing in some basis where the charged lepton mass matrix is not exactly diagonal 
j32| . This is a situation often encountered in realistic models f2]. Furthermore in GUT models 
it is often the case that, in the basis where the neutrino mixing is of the TB form, the charged 
lepton mixing matrix has a Cabibbo-like structure rather similar to the quark mixing and is 
dominated by a 1-2 mixing 0^ [3^, 



ei. 



CnE 


-Snse ^-^la 
"12 


o\ 


E 

E e^-^ia 


CgE 





12 


"^12 










1 J 



(4.1) 



where Cqe = cos , Sqe = sin 6*^ , and is a phase required for the diagonalisation of the 
charged lepton mass matrix [2]. The physical PMNS oscillation phase 6 turns out to be related 
to Af2 by 133], 

(4.2) 



Afs + vr. 



In this paper we assume that the neutrino mixing is accurately of the TB form. 



( 



\ 



'2 
3 
1 



1 

1 

V3 



\ 





1 

t 

V2 J 



p. 



Maj- 



(4.3) 



The physical mixing matrix is given by Eq j2.3| using Eq,4.3 and Eq j4.1[ The standard PDG 
form of the PMNS mixing matrix in Eq j2.4| requires real elements {UpMNs)n and {Upmns)i2 

. a.-\^ . o<.'2 

and this may be achieved by use of the phases in Puaj = diag{e'^~ , ,0). 
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It follows that {UpMNs)3ii {UpMNs)32 and {Upmns)33 are unaffected by the Cabibbo-like 
charged lepton corrections and are hence given by: 

\{UpMNs)3i\ = \iVlhi\ = ^, (4.4) 

|(^PAf7V5)32| = 1(^^321 = ^, (4-5) 
|(^PAf7V5)33| = 1(^^331 = ^- (4-6) 



The relations in Eqs,4.4 4.5 4.6 may be expressed in terms of the third family matrix element 
deviation parameters defined in Eq,2.8 as simply: 

C^ = 0. (4.7) 
Since these relations are all on the same footing, it is sufficient to discuss one of them only 



and in the following we choose to focus on Eq. 4.4 From Eq, 2.4 Eq. |4.4| may be expanded 
in terms of the standard mixing angles and phase leading to the so called mixing sum rules as 
follows: 

El = arcsin (V2 \s23S12 - si3C23Ci2e'^|') = 35.26°, (4.8) 



where we have assumed S23 = sin 6*23 = ■ This sum rule can be simplified further to leading 
order in S13, 

r2 = arcsin ( \/2 (S23S12 - S13C23C12 cos 6) ) « 35.26°. (4.9) 



From Eq,4.6 S23 = C23 = l/v2, hence to leading order, 

r3 = 9i2 - 013 cos{6) ^ 35.26°. (4.10) 

The last form of the sum rule was first presented in [2j , while all the forms can be found in |14] . 
We shall later study all three forms of the sum rules Fj, together with some related sum rules 
which we now discuss. 

Using the parametrization in Eq |2.5[ the sum rule in Eq,4.10 can be expressed in terms of 
the deviation parameters s, r and the Dirac CP phase ((5) [25], 



(Ti = rcos5 — s = 0. (4-11) 

To deal with issues of canonical normalisation corrections, the following sum rule has been 
proposed |19] . 

2 

(j2 = rcos5 + -a - s = Q. (4.12) 
3 

This sum rule was claimed to be stable under leading logarithmic third family RG corrections, 
although, as emphasized in [T9], it does not include the effect of the running of ^13 or r, whose 
inclusion introduces a Majorana phase dependence. Such effects will be studied numerically 
later. 



^This sum rule was derived from an expansion in 7712/7713, and the running of r was neglected because it is 
suppressed by an extra factor of -012 jmz compared to the running of s and a. 
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4.1.2 A GUT-Flavour Inspired Numerical Example 



In order to study the RG corrections and reliability of the various sum rules numerically it is 
necessary to define the GUT scale matrices rather specifically. In most of this paper we shall 
consider a numerical example inspired by the GUT-flavour models of [2} [6] , although in Section 
6 another numerical model will be considered leading to qualitatively similar results. Therefore 
in most of the remainder of this paper we shall take the right-handed neutrino Majorana mass 
matrix Mrji to be the diagonal matrix: 

/ 5.1 X 10-*^ \ 

Mrr = 7.05 X 10-9 Ms, 

V 1/ 

where M3 = W^^GeV. This is an example with light sequential dominance where the lightest 
right handed neutrino is dominant [18j. Ignoring RGE corrections to begin with, we find that 
precise tri-bimaximal neutrino mixing {612 = 35.26°, ^23 = 45.00°, = 0.00°) can be achieved 
with the Yukawa matrix: 

/ 1.0616676 0.001 \ 
Yi:R=\ e b (4.13) 

\ -0.9799e b C3 J 

where b = 8.125 x 10~^, e = 2.125 x 10"'^ and C3 = 0.5809. These parameters also lead 
to the following values for the neutrino masses: nii = 1.75 x 10~^eV^, 1712 = 8.67 x W~^eV, 
ms = 4.95 X IQ-'^eV, Am^^^ = 2.37 x IQ-^eV'^ and Am^^, = 7.52 x lO'^eV^. 

The low energy pole masses of the quarks are given as follows: = 1.22 MeV, rud = 
2.77 MeV, = 53 MeV, nic = 0.595 GeV, ruh = 2.75 GeV and mt = 163.6 GeV. In order to 
satisfy these values at low energy scale, REAP was used to perform the running of these masses 
from the Mz scale to the GUT scale and the resulting quark Yukawa matrices and Y^ at the 
GUT scale were taken as initial conditions for the running of the neutrino mixing parameters 
and sum rules from the GUT scale to Mz scale. 

The above parameter choice approximately satisfies the CSD conditions in Eq j3.7| However 
small corrections are used in order to achieve TBM neutrino mixing angles to 2 decimal places. 
If the CSD conditions were imposed exactly we would find instead 612 = 33.97°, $23 = 44.38°, 
^13 = 0.059° and 6 = 0° which are close to, but not accurately equal to, the TBM values. This 
is to be expected since the SD relations are only accurate to leading order in m2 / ms [18] . Since 
in this paper we are interested in studying the deviations from exact TB neutrino mixing due 
to charged lepton corrections and RG running we shall assume the matrices in Eq j4.13 rather 
than the CSD conditions as the starting point for our analysis. 

In order to study the effect of Cabibbo-like charged lepton corrections on the physical mixing 
angles where the neutrino mixing is precisely tri-bimaximal, we shall use the REAP package 
previously discussed. In order to use the REAP package it is convenient to work in the basis 
where the charged lepton Yukawa matrix is diagonal. Thus, assuming charged lepton corrections 
of the form of Eq j4.1[ the neutrino Yukawa matrix in the non-diagonal charged lepton basis must 
be transformed to the diagonal charged lepton basis according to: 

Y, ^ Y; = Ve^Y,. (4.14) 



Hence the original neutrino Yukawa matrix in Eq,4.13 must be rotated to the diagonal charged 



lepton basis according to Eq 4.14 
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Including the Cabibbo-like charged lepton corrections, physical tri-bimaximal mixing only 
holds when 0^ = 0. However according to the sum rules for Fj, certain combinatioms of mixing 
parameters sum to 35.262° for all values of the Cabibbo-like charged lepton corrections. This 
is illustrated in Tables J2] jS] where the values of the mixing angles together with the Dirac phase 
and the sum rules Ti, T3 at the GUT scale are presented for different values of 0^ and A^. 
Ti was found to be the most accurate sum rule at the GUT scale with a value of 35.262 ° exactly 
at all values of 6^2 and A{|. However the error in all the sum rules is less than about 0.1° in all 
the examples considered. 



^12 





1 


3 


5 


8 


012 


35.26 


34.648 


33.429 


32.216 


30.407 


013 


0.001 


0.708 


2.122 


3.534 


5.648 


023 


45.001 


44.997 


44.962 


44.892 


44.721 


d 





210.204 


210.82 


211.492 


212.672 


Ti 


35.262 


35.262 


35.262 


35.262 


35.262 


T2 


35.262 


35.26 


35.247 


35.217 


35.133 


Ts 


35.261 


35.26 


35.252 


35.23 


35.162 



Table 2: Values of the neutrino mixing angles 0i2, 0i3 and 023 together with 6 and the sum rules 
Ti, r2 and T3 at the GUT scale, at = 30° and tan{(3) = 50. All the angles are in degrees. 



^12 





7.5 


15 


30 


45 


012 


31.72 


31.752 


31.846 


32.216 


32.8 


013 


3.534 


3.534 


3.534 


3.534 


3.534 


023 


44.892 


44.892 


44.892 


44.892 


44.892 


5 


180 


187.9 


195.789 


211.492 


227.039 




35.262 


35.262 


35.262 


35.262 


35.262 


r2 


35.262 


35.259 


35.250 


35.217 


35.174 




35.254 


35.253 


35.248 


35.230 


35.208 



Table 3: Values of the parameters: 0i2, ^13,^231 ^ and the T sum rules at the GUT scale. These 
values are found in degrees at = 5 ° and tan(/3) = 50. 



4.2 More general charged lepton corrections including 

In the previous subsection we saw that the sum rules arising from Cabibbo-like charged lepton 
corrections are satisfied to excellent precision at the GUT scale, for the considered numerical 
example. In this section we introduce the case of non-Cabibbo-like charged lepton corrections. 
To be precise we shall consider more general charged lepton corrections given by. 



CgE 


— SnE e *'^12 





"12 


"12 




E e*^i2 


CnE 





12 


^12 










1 



/ 1 









CnE 

^Ke 

SgEC'^ 
"23 



23 





-SnE e~*'^23 
'^23 

CnE 
"23 



(4.15) 
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where we have now allowed both 0^ and to be non zero. The neutrino Yukawa matrix will 
be transformed to the diagonal charged lepton basis according to 



Y' 



Ve.Y,, 



(4.16) 



but now using the non- Cabibbo- like charged lepton rotations in Eq,4.15 After performing the 



charged lepton rotations in Eq,4.16, values for the mixing angles as well as the parameters 
given by Eq,2.8 can be calculated at the GUT scale. Of course in the present case of non- 
Cabibbo-like charged lepton corrections the third row deviation parameters ^i, ^2 and ^3 are 



all expected to be non-zero at the GUT scale. This implies that the sum rules given by Eq,4.7 
no longer apply in the case of charged lepton corrections with non-zero 0^ . The effects of non- 
Cabibbo-like charged lepton corrections on the deviation parameters E,i is displayed in Table |4] 
using the original neutrino Yukawa matrix as before, namely Eq 4.13[ but now with a small 
non-zero value of 0^ = 2°, and with different values of the new phase A^. 

Note that the effect of turning on the charged lepton correction 0^ will lead to a correction 



of the physical lepton mixing angle 622, but not 612 (to leading order) [2]. Therefore while the 
sum rules ri^2 and (T2 are violated by a non-zero the sum rules Ts and ai are both insensitive 



toOg 



■0 





161 


161 


161 





0.034 


0.034 


0.035 


30 


0.027 


0.031 


0.030 



Table 4: This table shows the values of |6M6I 
Cabibbo-like charged lepton corrections with 0^ : 
for different values of the phase A^. 



and 1^3 

= 5° 



^12 



at the GUT scale for case of non- 
= 30°, 6lf3 = 2° and tan(/3) = 50, 



5 Renormalization group running effects 

If TB neutrino mixing holds in the framework of some unified theory, then typically we expect 
Cabibbo-like charged lepton corrections leading to sum rule relations. However, as already 
indicated, such sum rules are only strictly valid at the GUT scale, and will be subject to RG 
corrections. In this section we now turn to a quantitative discussion of such RG corrections to 
the sum rules. For definiteness we shall assume the minimal supersymmetric standard model 
(MSSM), with a SUSY breaking scale of 1 TeV, below which the SM is valid. To study the 
running of the neutrino mixing angles and sum rules from the GUT scale to the electroweak 
scale, the Mathematica package REAP (Renormalization of Group Evolution of Angles and 
Phases) was used [T7] . This package numerically solves the RGEs of the quantities relevant for 
neutrino masses and mixing. It can be downloaded from http: / / www.ph.tum.de/ rge/REAP/. 
Mathematica 5.2 is required. 

^The insensitivity of the sum rule ai to 62Z is clearly seen numerically in Fig 6 (b). 
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5.1 Sum rules with Cabibbo-like charged lepton corrections 



5.1.1 Sum rules in terms of lepton mixing angles 

In this section, we study the RG running of the sum rules which result from Cabibbo-like charged 



corrections. The neutrino Yukawa matrix is taken to be of the form of Eq,4.13 as before. The 
RG change in the quantities, defined for a parameter P as A.P = Pmz ~ Pmgut^ calculated 
for the lepton mixing parameters and the Tj sum rules, and is presented in Tables, 5j6 From the 
results we see that the least precise sum rule F3 actually is subject to the smallest RG running 
since it does not involve ^23 which runs the most. 

The RG running of Fj is displayed in Fig{T]for tan(/3) = 50. The RG evolution of Fi and F3 
was also plotted at different values of tan(/3) as shown in Figj2j 
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+ 0.444 
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+ 0.151 
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+ 1.008 


+ 1.013 
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+ 7.453 


+ 2.126 


+ 1.181 


+ 0.62 


AFi 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.953 


AF2 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.954 


+ 0.958 


AF3 


+ 0.237 


+ 0.259 


+ 0.301 


+ 0.345 


+ 0.412 



Table 5: RG changes of the mixing parameters and sum rules Fi, F2 and F3 at 
tan{P) = 50. All values are in degrees 
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+ 0.444 


+ 0.432 


A013 
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- 0.091 


- 0.087 
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- 0.046 


A^23 


+ 1.009 


+ 1.009 


+ 1.009 


+ 1.008 


+ 1.006 


Ad 





+ 0.31 


+ 0.613 


+ 1.181 


+ 1.663 


AFi 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.953 


AF2 


+ 0.953 


+ 0.953 


+ 0.953 


+ 0.954 


+ 0.956 


AF3 


+ 0.362 


+ 0.36 


+ 0.357 


+ 0.345 


+ 0.326 



Table 6: RG changes of the neutrino mixing angles, the Dirac phase 5 and the sum rules Fi, 
F2 and F3 at 0^ = 5 ° and tan{f3) = 50. All values are in degrees. 



5.1.2 Sum rules in terms of TB deviation parameters 



In this subsection, for completeness we study the evolution of the TB deviation parameters 
Their RG evolution, for different values of 0^, is shown in Figlsl In FigJ4^we 



2.5 



defined in Eq 

display the evolution of the sum rules given by Eqs jlTTj [4T2| From Figjijit is seen that both 
(Ti, (J2 are precisely equal to zero at the GUT scale for = but differ by a tiny amount for 
612, ^12 7^ 0- 111 this numerical example it is apparent that the sum rule (J2 is slightly more stable 
than the original sum rule di, although there is not much more stability. This is a manifestation 
of the fact that (T2 does not take into account the running of r, which introduces an effect coming 
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Figure 1: Evolution of sum rules ri,r2, for Cabibbo-like charged lepton corrections for large 
tan(/3) = 50. Panel (a) shows the running at ^ = 5 ° and Aj^ 
evolution a t 



while panel (d) 



and 



0°. In panel (c) 



shows the evolution at 0^ 



5°, = 30°. Note how the graphs for Ti and r2 



0°, while panel (b) 
the running is at 



5 ° and Afa 



shows the 
15° 



completely overlap in (a) and (b) 
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Figure 2: Evolution of sum rules Ti, 
tan(/3). Running of Ti is shown in panel (a) 



for Cabibbo-like charged lepton corrections for various 
and Fa is shown in panel (b) at 0^ 



0° 



and different values of tan(/3). Note the expanded (and different) vertical scales used in these 
two figures; in all cases the corrections are less than one degree. 



from the Majorana phases which we have assumed to be zero in this example. Later on we shall 
discuss a numerical example with non-zero Majorana phases where the enhanced stability of C72 
will be more pronounced. 



5.2 Sum rules with more general charged lepton corrections including 9. 



E 
23 



Finally in this subsection we study the evolution of the parameters for the case of charged 
lepton corrections of the more general form in Eq,4.15 In Fig, 5 (a)| we show the RG running 
of the parameters ^i, ^2 and ^3, given in terms of the mixing angles in Eq. 2.9 for the case of 
Cabibbo-like charged lepton corrections. As expected, for Cabibbo-like charged lepton correc- 
tions, these parameters are exactly zero at the GUT scale for all values of 6^2 and A^, but then 



diverge from zero due to the RG corrections. In Fig. 5 (b) we now switch on the non-Cabibbo- 
like charged lepton corrections by a small amount corresponding to 0^ = 2 ° . In this case we 
see that the parameters ^1, ^2 and ^3 are all non zero at the GUT scale and deviate even more 
at low energies due to RG running. 

In FigjH] we show the running of the TB deviation parameters and the sum rules ai and cj2 
for the non-Cabibbo-like case with 0^ = 2 °. It is clear from panel (b) that the ai sum rule is 
still valid at the GUT scale even for a non-zero ^23, as remarked earlier. 
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Figure 3: Evolution of the deviation parameters r, s, a from the GUT scale to the electroweak 
scale without and with charged lepton corrections for large tan(/9) = 50. Panel (a) shows the 
running of the deviation parameters without any charged lepton corrections, 0^ = °, = °. 
Panel [(b)] shows the running of the deviation parameters in the presence of Cabibbo-like charged 
lepton corrections, 0^ = 5°, = 15°. 




Figure 4: Evolution of the sum rules cJi and ct2 from the GUT scale to the electroweak scale 
without and with charged lepton corrections for large tan(/3) 



without any charged lepton corrections, 0^ = 



and 



50. Panel 
° . In panel 



the presence of Cabibbo-like charged lepton corrections, 0^ = 5 ° and 
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Figure 5: The evolution of the third row deviation parameters ^i, ^2 and ^3 from the GUT 
scale to the electroweak scale at tan(/3) = 50 for Cabibbo-like and more general charged lepton 
corrections. Panel (a) shows the result with Cabibbo-like charged lepton corrections with 0^ = 

30 ° . Panel (b) is for the case of more general charged lepton 
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Figure 6: Running of the TB deviation parameters and their related sum rules from the GUT 
scale to the electroweak scale in the presence of more general charged lepton corrections, 0^ — 
5 °j ^12 = 30 °, Af3 



30°, 



23 



2 ° and tan(/3) = 50. Panel (a) 



shows the evolution of the TB 
deviation parameters. Panel |(b) | displays the sum rules ai and (T2. Note that cti = at the GUT 
scale even in the presence of the more general charged lepton corrections. 



16 



5.3 Sum rules with non-zero Majorana phases 



So far we have presented results for a particular example with zero Majorana phases. In this 
section, we present the running of the sum rules and the TB deviation parameters where the 



neutrino Yukawa matrix is taken to be similar to Eq,4.13 with the same values for b, e and C3 
but with non- zero Majorana phases: 





-1.012ee^^i 



0.972826e''^2 




(5.1) 



where we shall take the phases to be 61 = 120° and 62 = Q0°. The right-handed Majorana mass 
matrix is as before. The numerical value of the Yukawa couplings has been changed slightly to 
compensate for the non-zero phases in order to once again yield exact tri-bimaximal neutrino 
mixing at the GUT scale. 

In Fig. [7] we show results for the running of the sum rules and for the deviation parameters 
r, s, a for the above example with non-zero Majorana phases. In this example the (T2 sum rule 
is much more stable than cji as clearly shown in FiglT] (a) This shows that the question of the 



stability of the sum rule (T2 is dependent on the choice of Majorana phases via the running of 
r. In particular with this choice of Majorana phases the deviation parameters s, a and r all run 
less as shown in Figj7| [(b)| compared to the previous case with zero phases (Fig, 3 (a) ). 

The Fj and sum rules also change with the Majorana phases turned on but not as much 
as ai sum rules. For instance, at 0^ = 5° and = 0", we find that Fi and F2 get smaller by 
0.05 degrees at the Mz scale compared to the case where the phases are zero. F3 on the other 
hand gets larger by about 0.1 degrees. At = 5° and Afa = 30°, ^1 and 6 get smaller by 
about 0.001 to 0.003 compared to the zero phases case whereas ^3 gets larger by 0.006. 




Figure 7: Sum rules with Majorana phases. The running of cTj sum rules and TB deviation 
parameters are for Of^ = 0°, Xf^ = 0°, tan(/3) = 50, 61 = 120° and 62 = 60°. 
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6 Model dependence of the results: heavy sequential dominance 



So far all the numerical results have been based on a particular example inspired by the models 





of [21 [6], namely the case where the GUT scale neutrino Yukawa matrix has the form in Eq,4.13 
or the closely related form in Eq j5.1| with non-zero Majorana phases. In these examples the 
dominant contribution to atmospheric neutrino mass is coming from the lightest right-handed 
neutrino via the see-saw mechanism, a situation known as light sequential dominance (LSD) 
|18j . In order to test the generality of the results in this section we consider a quite different 
example in which the dominant contribution to the atmospheric neutrino mass is coming from the 
heaviest right-handed neutrino via the see-saw mechanism, a situation known as heavy sequential 
dominance (HSD) [18j. This example is chosen since it the most qualitatively different to the 
example of LSD considered previously, yet despite this we shall see that the numerical results 
for the corrections to TB mixing are qualitatively similar to those encountered previously. This 
gives us some confidence that our results and conclusions are not restricted to the particular 
numerical example studied but are in fact applicable to a large class of see-saw models based on 
hierarchical neutrino masses. 

In the HSD example considered here the right handed neutrino Majorana matrix as well as 
the neutrino Yukawa matrix are given by the following equations: 

-6 

Mrr 

where M3 = W^'^GeV. Ignoring RGE corrections to begin with, we find that precise tri- 
bimaximal neutrino mixing at the GUT scale {9'(2 = 35.26°, 9!^^ = 45.00°, e'(.^ = 0.00°) can be 
achieved with the Yukawa matrix: 

/ LOOl X lO""^ 1.0036 6 \ 
Y^j^ = b -1.0013 e (6.1) 

\ 2.992 X 10^-^ b e / 

where h = 2.401 x 10^^, e = 0.677. These parameters also lead to the following values for the 
neutrino masses: Amlt^ = 2.47 X W'^eV^ and Am^^^ = 7.53 x lO^^eF^ which are well within 
the allowed experimental ranges. 

Note that in the case of HSD the Yukawa couplings present in the neutrino Yukawa matrix 
are larger than the previous case especially e which we take to be 0.677 compared to 2.125 x 10~'* 
in the previous example. Furthermore there are similarly two large Yukawa couplings in the third 
column of the Yukawa matrix. Moreover the heaviest RH neutrino associated with these large 
Yukawa couplings has a mass well below the GUT scale leading larger threshold corrections 
coming from it. 

We assume charged lepton corrections of the form of Eq |4.1[ the neutrino Yukawa matrix in 
the non-diagonal charged lepton basis is then transformed to the diagonal charged lepton basis 
according to Eq |4.14 Using the REAP package, the running of Lj sum rules was performed 



from the GUT scale to low energy scale and the results are shown in FigurejS} As shown in 
this figure, despite the larger threshold corrections, for tan(/3) = 50, the RG running of is 
still small (about 0.4°) whereas that of Ti and r2 is about 1.3°, compared to the results shown 
in Figurej2] (nearly 1°). This suggests that, qualitatively, the results obtained for the previous 
numerical example inspired by the GUT models in [21 16] are expected to have wide applicability 
beyond the specific example considered. 
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Figure 8: Sum rules with heavy sequential dominance. This figure shows the evolution of the 
Ti sum rules for tan /? = 50 (panel (a) ), the running of Ti and Ts in terms of tan /3 (panels |(b)[ 



(c) ) from the GUT scale to the electroweak scale for the case of Cabibbo-like charged lepton 
corrections with flfa = 5 °, Afg = °. 
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7 Justification of the numerical approach 



The results in this paper have been based on a numerical evaluation of the RG corrections using 
the REAP package. The reasons why we have chosen to study these corrections numerically, 
rather than using analytic estimates [17] are as follows. 

The first reason we follow the numerical approach is that, as we show, some analytic esti- 
mates of RG effects which have ignored the effects of phases are unreliable. For example, the 
main purpose of this paper is to find out precisely how large the RG corrections are to sum 
rule relations which have been proposed in the literature (see section 4.1.1). Although the RG 
corrections to such sum rules are expected to be small, they are certainly not negligible com- 
pared to the expected precision of future neutrino experiments, and indeed this prompted the 
introduction of the modified sum rule in Eq j4.12 where the extra term compared to Eq,4.11 



was supposed to take into account the RG corrections [19]. However, it turns out that the new 
analytic term, which ignores the effects of phases, is too simplistic. Indeed the numerical results 



in Figures 4a and ^a^ clearly show that the extra term included in the analytic estimate of the 



RG correction in Eq, 4.12| does not capture the phase dependence of the RG correction to the 



original sum rule in Eq,4.11 The numerical study in this paper has highlighted the shortcom- 
ing of analytic estimates of the RG corrections to sum rules which do not include the phase 
dependence. 

The second reason we follow the numerical approach, rather than an analytic approach, is 
that for some of the cases studied the analytic approach is simply not applicable. The usual 
analytic approach is based on the assumption that only the third family charged lepton and 
neutrino Yukawa couplings are taken into account (while many analytic studies ignore neutrino 
Yukawa couplings and threshold effects altogether). Whilst the approximation of keeping only 
third family Yukawa couplings is sufficient for some models, for example the LSD class of models, 
it is certainly not sufficient for all classes of models. For example the HSD case that we also 
study involves two large neutrino Yukawa couplings, and the analytic estimates in [19j do not 
directly apply to this case. 

The third reason for following a numerical approach is a purely quantitative one, namely, 
even for the cases where the analytic approach is reliable and applicable (and we have already 
seen examples in the previous two paragraphs when it is neither) we would like to obtain the 
best possible estimate of the RG corrections which are the main focus of this paper. If the sum 
rules are to be confronted with experiment, it is important to have a reliable quantitative handle 
on the RG corrections, and for this purpose it is necessary to go beyond the leading log analytic 
approximation, as we now discuss. 

In order to investigate the quantitative accuracy of the analytic approach, in the remainder 
of this section we shall compare the analytic estimates of the RG effects for the LSD example 
studied earlier in the paper. For this purpose it is sufficient to switch off the charged lepton 
corrections and study the RG corrections to the neutrino mixing angles using the analytic 
approximations in [T7] which we then compare to the numerical results we obtained earlier in 
the paper, and which we also summarize here for convenience. In order to estimate the RG 
corrections to the mixing angles, following [TT] it is assumed that the (3,3) matrix elements 
govern both the charged Yukawa matrix (V^) and the neutrino Yukawa matrix (l^*^)!^ in the 
flavour basis in which the charged lepton mass matrix is diagonal. Taking ~ diag{0, 0, i/r) 
and Y'^ ~ diag(0, 0, Uus), one finds, to leading log approximation, that there is a single parameter 

^We have already noted that for some models such as HSD this is not the case for the neutrino Yukawa matrix. 
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which governs the RG corrections to all the mixing angles given by |19] : 



' 8^2^ Mz 87r2 



^-In 



Mgut 



(7.1) 



Assuming tri-bimaximal neutrino mixing at the GUT scale, the low energy scale parameters are 
then given approximately by: 



RG 



1 ry^RG 



^RG 
3 ma' 



(7.2) 



We now apply the above analytic formalism to the LSD model defined in section, 4. 1.2 



and subsequently studied numerically in this paper. In this model from Eq,4.13 we see that 
Vus = C3 = 0.58 at the GUT scale. We also find Ur = 0.33 and the mass ratio m2/m^ = 0.16 
for the case tan(/3) = 50. Using these values, the mixing angles and the quantities (Tj) can be 
estimated as shown in Table jT] where the analytic estimates are shown alongside the numerical 
results for comparison. 



Parameter 


^12 


^23 


^13 


Ti 


r2 




Analytic (°) 


35.59 


45.69 


0.15 


35.94 


35.94 


35.44 


Numerical (°) 


35.65 


46.00 


0.15 


36.21 


36.21 


35.49 



Table 7: This table shows a comparison between the analytic and numerical results for the RG 
corrections to neutrino mixing angles at the Mz scale, assuming that they take the precise TB 
mixing values at the GUT scale, for the LSD model described in the body of the paper with 
tan(/3) = 50. 

It is interesting to compare the analytic results to the numerical results in Table [7] for the 
neutrino mixing angles at the Mz scale, assuming that they take the precise TB mixing values 
at the GUT scale and setting all charged lepton corrections to zero, for the LSD model described 
above. The results show that the numerical estimate of O^^ (which is equal to zero at the GUT 
scale) is very accurately reproduced by the analytic approximation (indeed there is no difference 
to 2 d.p.), and the RG correction to 9 12 is also well reproduced with the analytic estimate 
underestimating the correction by only 0.06 degrees. However the results also show that there 
is a significant underestimate of ^23 with the analytically estimated value at the Mz scale being 
less than the numerical value by about 0.3 degrees, resulting in the analytically estimated values 
for Ti and r2 being less than the numerical values by about the same amount (0.3°). From the 
point of view of the effects studied in this paper (for example note the precision of the scales 



shown in the results in Figure la) an error of 0.3 degrees is undesirable and we would not wish 
to compromise the results in this paper by being subject to such unnecessary errors incurred by 
the analytic approach. 

Finally we remark that the origin of the discrepancy between the analytic estimates and the 
numerical results, for the cases where the analytic approach is reliable and applicable, is due 
to the fact that the analytic estimates are based on the assumption that the Yukawa couplings 
Ur and Ui,^ are fixed at their GUT scale values and do not run, whereas the numerical results 
allow for the co-running of all the Yukawa couplings in the matrix (including the second family 
Yukawa couplings), with the leading logs being effectively re-summed. 
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8 Summary and Conclusions 



In this paper we have analyzed the effects of charged lepton corrections and RG running on 
the low energy predictions of theories which accurately predict tri-bimaximal neutrino mixing 
at the high energy scale. In GUT motivated examples the charged lepton corrections are often 
Cabibbo-like and in this case the effect of charged lepton corrections leads to a range of neutrino 
mixing sum rules at the GUT scale, given by the Tj sum rules in Eqs, 4.8|4.9|4.10 as well as 



the (Tj sum rules expressed in terms of the deviation parameters in Eqs,4.11 and 4.12 We 



have studied the RG running of such sum rules numerically for a specific numerical example 
inspired by the GUT models in [21 [6], corresponding closely to CSD with LSD. Our results 
indicate small but measurable effects for the two examples studied. For example the sum 
rule arising from Cabibbo-like charged lepton corrections (due to 0^ corrections) which at the 
GUT scale corresponds to 612 — ^13 cos((5) 35.3° becomes renormalized by about 0.4° even for 
large tan P = 50. We have also considered the effect on charged lepton corrections coming from 
non- Cabibbo-like charged lepton corrections (due to non-zero and confirmed that the sum 
rule (Ti is insensitive to 

Even for a particular class of numerical model, such as the GUT-flavour inspired LSD model 
considered, the numerical results will depend in general on the choice of Majorana phases for 
that model. We have seen that switching on these Majorana phases can alter significantly the 
running of the TB mixing deviation parameters r,s,a as well as the sum rules such as cij. For 
example the sum rule (T2 which includes the leading logarithmic RG corrections due to the 
running of s and a, will have a Majorana phase dependence via the running of r which was 
neglected in the derivation of (T2 [1.9,. Thus, the relative stability of C72 as compared to ai turns 
out to be a Majorana phase dependent question. 

Although most of the numerical results are based on a particular GUT-flavour motivated 
LSD type of model, we have also considered similar results for a completely different type of 
model based on HSD. Overall we have found that the RG running effects are quite small in 
both cases which suggests that qualitatively similar results will apply to other models based 
on the Minimal Supersymmetric Standard Model, extended to include the see-saw mechanism, 
with hierarchical neutrino masses. However, we repeat our caveat in the Introduction that in 
this paper we have only considered deviations from TB lepton mixing due to the combination 
of charged lepton corrections and RG running, and that in general there will be other sources of 
deviations that we do not consider. For example, exact CSD [2] itself will lead to some deviations 
since it does not predict precisely TB neutrino mixing due to corrections of order 1712 / fTE\ . 
Another example of deviations to TB mixing are the canonical normalization effects discussed 
in [IH]. In general corrections to the TB mixing in the neutrino sector, in the flavour basis, are 
highly model dependent and have been studied phenomenologically in |20j . 

Finally we note that all the results in this paper are based on a numerical approach rather 
than an analytic approach. We have given three reasons to justify the use of the numerical 
approach: flrstly that analytic results which ignore the effects of phases are not reliable (for 
example the C72 sum rule); secondly that analytic results may not be applicable for all cases 
of interest (for example HSD with two large neutrino Yukawa couplings); and thirdly that the 
quantitative precision of the analytic approach is not in general sufficient for our purpose here (in 
particular we have shown that the RG correction to 823 may be signiflcantly underestimated in 
the analytic approach). Therefore throughout this paper we have adopted a numerical approach, 
exploiting the convenient REAP package. 
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To conclude, we have assumed precise TB neutrino mixing at the GUT scale, and investi- 
gated the deviations due to both charged lepton corrections and RG running. In GUT-flavour 
inspired models the charged lepton corrections are expected to resemble those of the quark mix- 
ing matrix, and lead to neutrino mixing sum rules which are valid at high energy. We have 
studied numerically a variety of such neutrino mixing sum rules and shown that, in the types 
of GUT-flavour models in which the sum rules emerge in the first place, they are subject to 
only mild RG corrections (less than one degree for the cases studied) in evolving them down to 
low energy. However even such small corrections as described in this paper will nevertheless be 
important when comparing the neutrino mixing sum rules to the results of future high precision 
neutrino oscillation experiments |34j . 
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